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We analyse the three-particle scattering continuum in 
quasi one dimensional integer spin Heisenberg antiferromag- 
nets within a low-energy effective field theory framework. We 
exactly determine the zero temperature dynamical structure 
factor in the 0(3) nonlinear sigma model and in Tsvelik's 
Majorana fermion theory. We study the effects of interchain 
coupling in a Random Phase Approximation. We discuss the 
application of our results to recent neutron-scattering exper- 
iments on the spin-1 Haldane-gap material CsNiCla. 
PACS: 75.10. Jm, 75.45. +j, 75.50.-y 



I. INTRODUCTION 

In recent inelastic neutron scattering experiments |^] 
on the quasi one dimensional spin-1 Heisenberg mag- 
net CsNiCla the existence of incoherent multiparti- 
cle scattering continua in the one dimensional phase was 
investigated. If was found, that close to the antiferro- 
magnetic wave vector (along the chain direction), there 
is significant spectral weight above the coherent magnon 
peak. 

Motivated by these experimental results we determine 
the dynamical structure factor for weakly coupled inte- 
ger spin Heisenberg chains in a low-energy effective field 
theory framework. In particular, we calculate the ratio 
of spectral weights of multiparticle scattering continua to 
the coherent magnon peak. 

An appropriate model Hamiltonian for CsNiCls is H 



H 



(1) 



<ij>' 



where the first and second sums are over nearest neigh- 
bour spins along and between the chains respectively. 
The exchange constants are estimated to be J « 2.8meV, 
J' « 0.045meV and the single-ion anisotropy is estimated 
to be Z) « — 0.004meV Above the ordering tempera- 
ture of about 4.84K, CsNiCls is considered to be a good 
realisation of a one dimensional Haldane-gap system, al- 
though there is a sizeable dispersion perpendicular to the 
chain direction due to the interchain coupling. As D is 
very small, we neglect it from now on. 

As a first approximation we can neglect the effects due 
to J', so that we arrive at the purely one dimensional 
Heisenberg Hamiltonian 



H 



ID 



n+1 



(2) 



The dynamical susceptibility for (g) has been calculated 
numerically for a chain of 20 sites by Takahashi ||^. He 
found, that about three percent of the total intensity 
above the antiferromagnetic wave vector are due to in- 
coherent multiparticle scattering continua. A recent dy- 
namical density matrix renormalisation group calculation 
an large systems of up to 320 sites quotes a result of 
about 2.5 percent for this quantity. The purpose of 
the present work is to determine the incoherent contri- 
bution to the structure factor by analytical means in the 
thermodynamic limit. At present, the best way to do this 
is by using a low-energy effective field theory description 
of (^ and taking the interchain coupling in (|^) into ac- 
count perturbativcly. The dynamical magnetic suscepti- 
bility for a spin chain is given by 



poo pc 

X"^{cu,q) = — / dx / 



dt e"^*-*«^ 



X ([5"(t,a;), 5^(0,0)]). 



(3) 



Note that we use units in which h — I. The dynamical 
structure factor is obtained as (see e.g. |5|) 



(4) 



Due to the spin-rotational symmetry of (|^), we can re- 
strict our attention to the case a — /3 — z. We denote 
the corresponding susceptibility and structure factor by 
x(ti',q) and S{uj,q), respectively. Below we will calculate 
the structure factor in two different low-energy effective 
field theories. We concentrate on the region g ~ tt, which 
is of experimental relevance. It is known that in this re- 
gion there exists a coherent one magnon excitation with 
a gap A and incoherent three, five, seven etc magnon 
scattering continua. We constrain our analysis to the 
case of three magnons, which gives the dominant contri- 
bution. We also investigate the effects of the coupling 
between chains, but do not take into account the single- 
ion anisotropy. We note that the analysis of section |^ 
can be easily extended to the case D ^ 0. 

The outline of this paper is as follows. In section || 
we determine the dynamical structure factor of the one 
dimensional model (^ for general integer spin S in the 
framework of the O (3) nonlinear u- model (NLcrM) de- 
scription. In section |l| we study the effects of the cou- 
pling between chains in (|^), using the NLctM results as 
an input. In section |^ we determine the structure factor 
within a second low-energy effective field theory descrip- 
tion, which holds for 5 = 1 in the presence of a (strong) 



1 



biquadratic exchange interaction Q between spins in ad- 
dition to (0). Although the structure of the excitation 
spectrum as a function of wave number along the chain 
direction is essentially the same as in the NLcrM , the 
structure factor turns out to be rather different. In sec- 
tion we discuss temperature effects and in section VI 
we summarise and discuss our results. 



II. 0(3) NONLINEAR SIGMA MODEL 

At energies much smaller than the exchange J, the 
lattice model (||) (for integer spin S) can be approximated 
by a field theory, the 0(3) NLcrM (see e.g. [|,|). The 
Lagrangian density of the NLcrM is given by 



dr 



dt 



dr 



dx 



fh ■ rfi = \ , (5) 



where m is a three-component vector and v is the magnon 
velocity. The relation between lattice and field theory 
variables is given by |f7|,p| 



S(x) w Si-lf'^^mlx) + —m{x) x 

vg at 



(6) 



where ag is the lattice spacing. In order for the NLcrM 
to give an accurate description of the low-energy physics 
of the spin chain, the spin S is supposed to be large. 
However, previous investigations indicate that even 
the case = 1 is described rather well by the NLcrM. 
Using (^), it is possible to determine dynamical correla- 
tion functions within the framework of the NLcrM . The 
q « behaviour of the dynamical structure factor was 
determined in this way by Affleck and Weston some time 
ago jf^] . For g w the NLcrM predicts an incoherent two 
magnon scattering continuum. As the scattering inten- 
sity is proportional to g^, this continuum is very difflcult 
to observe experimentally. 

Here we are interested in wave-number transfers close 
to the antiferromagnetic point g « Using (^, we see 
that the dynamical susceptibility is given by 



_„C'2 rca poo 

X(cj,<?) = / dx dte'"*-*9" 

ao 7_oo Jo 

x([m3(t,a;),m3(0,0)]) 



(7) 



The two-point function {ma{t,x) TOa(0, 0))has been cal- 
culated using the formfactor approach [l|] in |]l2| , |l3| . 
Let us briefly review some relevant formulas. The exact 
spectrum of the 0(3) NLcrM consists of three massive 
magnons Aa, a = 1,2,3, that form the vector represen- 
tation of 0(3). We denote the magnon mass gap by M. 
Due to factorisability of the exact scattering matrix , 
multi magnon scattering states form a basis of the Hilbert 
space. Let us introduce some notations. We parametrise 
energy and momentum of one magnon states in terms of 
a rapidity variable 9 



El = Ai'w^coshe , Pi = Mvainht 



(8) 



A scattering state of N magnons with rapidities {9j} and 
0(3) indices {sj} {sj — 1,2,3) is denoted by 



1^*1, 6*2, ■ ■ • , 6'Ar)ei,e2,...£N ■ 

Its energy and momentum are 

N N 

En = Mv^ ^ cosh 9^ , Pn = Mv ^ sinh Oj 
The resolution of the identity is given by 



(9) 



(10) 
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n=0 {c,} 



'3 = 1 



d9^ 
An 



X 1^1, . . . , 9n)ei,...e„ (6'i, . . . , 0„| 



(11) 



The two point function of some operator O can now be 
expressed in the spectral representation as 



1 rift 

{O\t,x)O{0,0))=j:^J U^l^o'iO^ 

n=0 '•'~°°j=l 
/ n n 

X exp I —it v'^M cosh 9^ + ix vAI sinh 6 
Here 



(12) 



fc=i 



k=l 



\:f^"\Oi, 



^|(0|Ot(0, 0)1^1,..., Oei...eJ 



(13) 



The matrix elements for several operators O in (|3|) have 
been determined in [|l^,^. In order to calculate the 
structure factor, we are interested in matrix elements of 
(any of the components of) the fundamental field of the 
NLcrM , e.g. O = m^. The summed absolute values of 
the first few matrix elements for this case are 

|.F(2^-)(0)P=O, fc=l,2,..., 
|.F(i)(0)|2 = Z , 

\T^'\9i,92,9s)\' = Z7r' 

n 



12Tr2 + 2 ^(6lj-6lj 



fc<; 



{Bu-9iY{\9u-9if + \n^) 



'j<i 



tanh ■ 



(14) 



Here the overall factor Z is due to the field renormalisa- 
tion in the NLcrM . In order to determine the dynamical 
susceptibility, we now use (|l|) and (H) in (0) and then 
carry out the x and t integrations. In this way we arrive 
at 
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f{y + z)f{y~z)f{2z) 



127r2 + 4(3z2 + y^)] 



^2 _ ^2g2 _ /^2 _ 4^2 coshz(cosh2 + cosh?/) + is 
' contributions from 5, 7,9... particles, (15) 



where Z' = S'^Zv/aa, A = Mv'^ and 



fiz) = [tanh(z/2)]^ 



z2 + ^2 
z2[z2 + 47r2] 



(16) 



The structure factor close to the antiferromagnetic wave 
number is obtained from this using (^) 



S(uj,q) 



Z' 



S{lu - ^„2q2 + A2) 



2VwV + A2 
-^Jdz [3^2 + 3^2 ^ y2] /(2;^)/(^ + r)/(^ _ Y) 

1 

[(u;2 - w2g2 _ A2 - 4A2 cosh^ z)2 - 16A4 cosh^ z] 
contributions from 5, 7,9... particles. (17) 



arccosh^Y^ and Y 



arccosh 



2;"'-l-4cosh"' 



Here zo - ......... ^ . - ......... ^^^^^^^ 

where — (uj^ — v'^q^A'^. An important differ- 
ence between (p^ ) and (|l5| ) is the way, in which states 
with 3,5,7... magnons contribute. It can be easily seen 
from the upper limit of integration, that three magnon 
states contribute to the structure factor only if s2 
ti;2 _ ^,2^2 ^ 9A2, i.e. above the three magnon thresh- 
old. Similarly, 2n + 1 magnon states only contribute 
if s > (2n -f 1)A. In other words, ( p7|) is exact as 
long as s < 5A. The situation is quite different for 
Rex(w,'Z): here multi magnon states contribute for all 
values of s. However, their contribution is negligible at 
small s (see e.g. H^). The first contribution in (17) is 
due to the coherent one magnon states with dispersion 
uj — -y/ u2g2 _|_ ^2^ which is known to be a good approx- 
mation to the lattice-dispersion as long as q is sufficiently 
small (recall that in our notations q denotes the deviation 
from Tr/ao)- The total spectral weight cannot be calcu- 
lated within the NLctM framework, so that one has to 
resort to a direct numerical analysis of the quantum spin 
chain. For the spin-1 case this yields Z' « 1.28^, 
V ~ 2.5JaQ. 

The remaining integral in the three magnon contribu- 
tion to ( p^ ) can be evaluated numerically. The result is 
shown in Fig . |[ 
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FIG. 1. NLaM result for the three magnon contribution to 
the dynamical structure factor close to the antiferromagnetic 
point. Her e M is the magnon mass-gap, v the spin velocity 
and s — yVg^M-A^. 

The behaviour of the structure factor just above the 
three magnon threshold can easily be determined by 
Taylor-expanding the integrand in (|lj) . The result is 



S{uj,q) 



Z' 



3A]^ 



0< — -3<1 . (18) 



The ratio of total spectral weights of the one magnon (/i) 
and three magnon contributions can be calculated 
as well. Numerically we find that that ratio of spectral 
weights at wave number tt/co is roughly equal to 

^3(^) 



0.02 . 



(19) 



We note that a sizeable fraction of the three magnon 
spectral weight is located at very high energies, where the 
field theory does not relate to the lattice spin model. If 
we restrict uj to be smaller than twenty times the magnon 
gap, the ratio ( [l6| ) diminishes to about 0.013. 

The NLctM result is somewhat at odds with the nu- 
merical results on the spin 1 chain. The most likely 
reason is that the NLcM does not work as well for 5=1 
as previosuly thought. Irrespective of the one percent 
difference between NLcrM and numerical results it seems 
clear that the three magnon continuum of a single spin-1 
chain is extremely weak! In quasi one dimensional mate- 
rials like CsNiCls such a faint contribution would hardly 
be measurable experimentally as it would go under in 
the inevitable errors due to e.g. background subtraction. 
Another relevant quantity as far as the experiments of 
||l| are concerned, is the ratio of three and one magnon 
spectral weights for q ^ 0. We find that e.g. 

/3([l±0.2]7r) 



0.04 



(20) 



/i([l±0.2]7r) 

The increase in ( |20| ) as compared to ( |l9| ) is mainly due 
to the decrease in spectral weight of the single magnon 
peak 
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/i([l±0.2]7r) 



0.45 



(21) 



Let us now determine the real part of the dynamical 
susceptibility, which we will need in section III. Neglect- 



ing contributions of more than three magnons we arrive 
at the results shown in Fig. 0. 
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FIG. 2. Minus the real part of the dynamical susceptibility 
calculated in the framework of the NLcrM . Contributions of 
states with more than three magnons are neglected. The inset 
shows the contribution due to three magnon states. 

We note that the only singularity in the real part is at 
s = A. 



III. COUPLING BETWEEN THE CHAINS 

Materials like CsNiCla are only approximately one di- 
mensional. There always is a small coupling between the 
individual chains, which we will now take into account 
within the framework of a Random-Phase Approximation 
(RPA) . Denoting the Fourier transform of the interchain 
coupling matrix elements by J(fc), the three-dimensional 
dynamical susceptibility in RPA is simply given by 1 15 



(22) 



Note that this generally is a matrix equation. In our case 
( p^ ) reduces to a scalar equation due to SU(2) symmetry 
of the full, three-dimensional Hamiltonian. We now can 
insert the results ( p^ ) and ([iT^ fo r the dynamical suscep- 
tibility of a single chain into (|22|) and obtain in this way 
the dynamical structure factor of weakly coupled chains 
for (7 w — 



S{uj,q, k) = 



S{uj,q) 



1- J(fc)x(^,g)P 



(23) 



Note that ( p3| ) involves both the real and imaginary parts 
of the one dimensional dynamical susceptibility. The 



magnon dispersion close to the antiferromagnetic point 
along the chains is easily extracted from the poles of (|2^ ) 



= v^q^ + A^ + Z'j{k) 



(24) 



The coherent one magnon contribution to the structure 
factor is found to be 



S{uj,q,k)\^ 



Z' 



magnon 



2Jv^q^ + + Z'J{k) 



X SiuJ- \^v^q^ + A2 + Z'J{k) 



(25) 



We see that the total spectral weight due to one magnon 
processes depends on the transverse wave-number trans- 
fer k. 

From now on we consider the particular case of wave- 
number transfer k — {rj, rj, 0) along the (1,1,0) direction in 
CsNiCla, which is of direct experimental relevance [|l|,||. 
The rj dependence of J{k) is of the form 



J(fc) = 2 J' (cos47r?7 + 2cos27r?7) 



(26) 



We note that at the special point ko — (0.19,0.19,0) we 
have J{k) = 0, so that (in our approximation) we are 
dealing with an ensemble of uncoupled chains. By fitting 
( p^ ) with ( p6|) to the experimentally observed magnon 
dispersion along the (1,1,0) direction pi, we obtain 



1.32meV 



Z'J' w 0.475meV2 



(27) 



The resulting value of J' w 0.05mcV is by construction 
close to The variation of the total spectral weight 
/i(7r, 77) of the one magnon peak with 77 (|25| ) (with fixed 
wave number n/ao along the chain) is shown in Fig.^ 




0.40 



FIG. 3. Ratio of the total spectral weights of the one 
magnon peak Ii{n,r]) for wave- number transfer rj along the 
(1,1,0) direction to the "one dimensional" result 7i(7r,0.19). 
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The three magnon contribution to the structure factor 
is of the form 



S{uj,q,k)\ 



S{uj,q)\ 



3 magnons 



3 magnons 



\l~Jik)x{u;,qW 



(28) 



Here we used the fact that the real part of the one di- 
mensional dynamical susceptibility did not exhibit any 
singularities. A result of this is that the three magnon 
continuum for fixed rj starts at three times the magnon 
gap only at the special point 77 = 0.19. 

From the analysis in section |l| we know that the dom- 
inant contribution to x{^^ q) for to > 3A is due to the 
real part of the one magnon contribution, which gives 
us an easy way to estimate the effects of the interchain 
coupling on the three magnon continuum 



S(uj,7r/ao,k)\„ « Sfw, 7r/ao)L 

^ ' ' ^ ' ^ \ 6 magnons ^ ' ' ^ ^ I J magnons 



1 - 0.547 



cos Anr] + 2 cos 2Trri 



1 



(29) 



In order to estimate the ratio of spectral weights of the 
three and one magnon states as a function of 77 we define 



20M 



duj S{uj, tt/co, k)\ 



3 magnons 



(30) 



Here the cutoff at twenty times the magnon gap has been 
chosen arbitrarily. Note that the introduction of a cut- 
off is necessary to ensure the applicability of the field- 
theory description. The ratio of the three magnon spec- 
tral weight I3 (tt, rj) to the one magnon spectral weight 
/i(7r, 77) for wave-number transfer rj along the (1, 1,0) is 
shown in Fig.^. 



0.025 



0.020 



" 0.015 



0.010 



0.005 



0.00 



0.10 



0.20 



0.30 



FIG. 4. Ratio of the total spectral weights of the three 
magnon continuum Is'{n,Tj) for wave-number transfer rj along 
the (1, 1,0) direction to the total one magnon spectral weight 
/i(7r,77). 



We see that in all cases only a very small fraction of the 
total spectral weight sits in the incoherent three magnon 
scattering continuum. This fraction is at most two per- 
cent, which is much smaller than what has been observed 
experimentally in CsNiCls [|| . We conclude that the cou- 
pling between chains cannot account for the observed in- 
tensity in the incoherent scattering continuum at n/ao 
as long as we use a NLcrM description for a single chain. 



IV. TSVELIK'S MAJORANA FERMION THEORY 

Apart from the NLtrM there is a second low-energy 
effective field theory for Haldane-gap systems, due to A. 
M. Tsvelik m&. It applies to spin-1 models of the type 



Sn^l — b I Sn ■ S, 



+ DiS^J^, (31) 



where « 1. We note that one of the effects of 7^ is 
to reduce the magnitude of the Haldane gap, which rules 
out a very large value of b in CsNiCls. Using the Bethe 
Ansatz solution of ( |3l] ) at the Armenian point b = 1, 
D — |]l^ , it is possible to derive jl^ the following low- 
energy effective field theory of three interacting Majorana 
fermions valid for small 1 — 6 and D 



maXaXa + 9aJnJn 



(32) 



Here a = 1,2,3, Xa are two-component Majorana (real) 
fermions, Xa = xllo, = ^"'^''XbluXc are the compo- 
nents of SU(2) currents and our conventions for gamma 
matrices are 



70 



71 



(33) 



In the absence of a single- ion anisotropy {D = 0), the 
model is SU(2) symmetric and all Majorana masses and 
couplings ga are equal. We will constrain our discussion 
to this case only. Following |l^ we will furthermore ne- 
glect the current-current interaction i.e. set ga = 0. This 
is certainly justified as long as the mass gaps are not too 
small and we are only interested in single-particle prop- 
erties as can be seen from a standard one-loop renormal- 
isation group calculation. However, it is presently not 
clear whether this remains a good approximation for the 
calculation of the three magnon scattering continuum. In 
fact, it was recently argued jisj that the current-current 
interaction may become rather important for the calcu- 
lation of finite-temperature properties of spin-1/2 ladder 
models, which have a very similar field-theory descrip- 
tion After setting = 0, (^) reduces to a theory 
of three noninteracting Ising models and exact informa- 
tion on correlation functions is available. The staggered 
components of the lattice spin operators are expressed in 
terms of order (a) and disorder (/i) operators of the three 
Ising models as 

(-l)^/"«S'^(t,a;) = VAa\t,x)fi^it,x)fi^t,x) , (34) 
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where the coefficient A is presently not known. The dy- 
namical susceptibility close to the antiferromagnetic wave 
number tt is thus given by 



x([(aVV')(i,2;),(aVV)(0,0)]) . (35) 



In order to determine (^5|) we therefore need to calculate 
{G^{t, x))'^Ga{t, x), where 



G^{t,x) = (a{t,x) ct(0,0)) , 
G^{t,x) = {fi{t,x) pi{0,0)) , 



are correlation functions of order and disorder operators 
in the two-dimensional Ising model. These are known ex- 
actly po| . For our purposes it is most convenient to work 
in the Minkowski-space spectral representation, which 
reads H 



°° 1 '"_ rid 



k=l 



k=l 



X exp 1 -it ^ v^M cosh 0k + ix'^ vM sinh 9k j , (37) 
where 



2n+l 

Wi) = i ' 



2n+l 



JJ^ tanh ■ 



t<j=i 



2n 



J^tanh^. (38) 



l<j=l 



After performing the Fourier integrals we arrive at the 
following representation for the dynamical structure fac- 
tor 



S{uj,q) 



vA 



+ 



4vA /"^f . tanh^ z + ^ [tanh z tanh tanh ' 



7r2A2 



dz- 



'0 ^(x2-l-4cosh2z)2-16cosh^z 
+ contributions from 5,7,9... particles , (39) 

where x^ = {liP — v'^q^)/A'^ and zq and Y are defined 
above. The one and three particle contributions in ( |39| ) 
give the exact result as long as s < 5 A. For frequen- 
cies just above the three-particle threshold the integral 
in (^9|) can be easily evaluated by Taylor-expanding the 
integrand 



S{uj,q) = 



Av UJ^ ~ v^q^ - 9A2 



(40) 



Eqn. ( ^0| ) agrees with a result obtained for correlation 
functions of the spin- 1/2 ladder [|l^,p2|. 



Performing the remaining integral in ( ^9| ) numerically, 
we obtain the three magnon contribution to (|3^) as 
shown in Fig. |^. 
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FIG. 5. Result for the three magnon contribution to the dy- 
namical structure factor close to the antiferromagnetic point 
in Tsvelik's Majorana fermion theory. H ere M is th e magnon 
mass-gap, v the spin velocity and s = v'^q^ -\- IS? . 

The ratio of spectral weights of three-particle and one- 
particle contributions to the structure factor ( p9[ ) for fre- 
quencies restricted to w < 20M is roughly equal to 



^i(^) 



0.17 . 



(41) 



Thus the three-particle scattering continuum is much 
stronger than in the NLcrM ! The coupling between chains 
can be taken into account in complete analogy with the 
NLcrM case. Furthermore it is in principle possible to 
analytically study the effects of the single-ion anisotropy 
in the Majorana fermion theory. 



V. TEMPERATURE 

An additional complication in CsNiCla is that the or- 
dering temperature 4.84K is of the same order of mag- 
nitude as the Haldane gap A « 15.4K. Clearly tem- 
perature effects will therefore not be negligible in the 
temperature range in which Neutron scattering experi- 
ments are conducted (5K-10K). One way of taking them 
into account would be to use recent results for finite- 
temperature correlation functions of the Ising model p^ ] 
in order to extend the analysis of section IV. This is 
nontrivial, so that we constrain ourselves to a rather pre- 
liminary discussion based on recent results of |]l8| . There 
it was shown, that the finite temperature lineshape of 



the one magnon peak at very low temperatures T <C A 
is given by 

S'(w,g) = yi'rcxp(-A/T) 
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(to - v/A2+t;2<j2)2 ^ (i.2rexp(-A/T))2 



(42) 



In other words, at temperatures much smaUer than the 
Haldane gap the Hneshape is Lorentzian. Although ( ^ ) 
does not really apply to CsNiCla because T = 0{A), it 
can be taken as an indication that in the experimentally 
relevant temperature range, a sizeable fraction of what 
used to be the one magnon spectral weight at T = may 
get transferred to energies significantly above the single 
magnon gap A. 



VI. SUMMARY AND CONCLUSIONS 

Using a low-energy effective field theory description 
(that needs to be supplemented by numerical results in 
order to fix overall normalisations) we have studied the 
dynamical structure factor for weakly coupled integer 
spin Heisenberg chains. The three particle incoherent 
scattering continuum calculated in the 0(3) NLctM is 
found to be too small to be observed in neutron scatter- 
ing experiments. 

On the other hand, Tsvelik's Majorana fermion theory 
yields a very strong three particle scattering continuum, 
which would be easily observable in experiments. The 
main difference between the NLtrM and the Majorana 
fermion theory is that in the latter the spin Hamiltonian 
contains strong biquadratic interactions of the form 



(43) 



The presence of such interactions in e.g. CsNiCls cannot 
a priori be ruled out. As a matter of fact we believe that 
spin-1 materials generically have at least small interac- 
tion terms of the form (^) . 

In the neutron scattering experiments a significant 
incoherent scattering continuum has been observed for 
several different temperatures. These findings are incom- 
patible with the NLcrM results. However, the results of 
section IV indicate that biquadratic exchange interac- 
tions lead to a transfer of spectral weight from the coher- 
ent magnon peak to the incoherent scattering continuum. 
One possible scenario for reconciling (field) theory with 
the experimental findings is thus to postulate the exis- 
tence of a sizeable biquadratic exchange interaction in 
CsNiCk. To the best of our knowledge no detailed (be- 
yond |2J]) theoretical study of the effects of biquadratic 
exchange interactions on the dynamical structure factor 
has been carried out. We believe that it would be inter- 
esting to do so. 

As far as the present analysis is concerned, it has sev- 
eral shortcomings that ought to be improved upon in 
future work. Firstly, our analysis is restricted to zero 
temperature. As we argued above, temperature effects 
are important in the experimentally realised range of pa- 
rameters. However, we believe it is unlikely that the find- 
ing of 0] are due to temperature effects only. Secondly, 



the analysis in the framework of the Majorana fermion 
theory did not take into account the current-current in- 
teraction. The results of section [V therefore should be 



regarded with some caution. It would be very interesting 
to carry out a systematic analysis of the effects of the 
current-current interaction. 

Finally we would like to mention that the analysis of 



section IV applies with minor modifications to the case of 
the spin-1/2 ladder as well. This is because the effective 
field theory description is very similar |19|| . 



While this paper was being written, a preprint (cond 
mat/990743ll ) by M.D.P. Horton and I. Affleck appeared 



in which the three magnon contribution to the dynami- 
cal structure factor in the NLctM was calculated. Their 
results have a strong overlap with our section ||. 
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